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Motivated by experiments in which single-stranded DNA with a short hairpin loop at one end 
undergoes unforced diffusion through a narrow pore, we study the first passage times for a particle, 
executing one-dimensional brownian motion in an asymmetric sawtooth potential, to exit one of the 
boundaries. We consider the first passage times for the case of classical diffusion, characterized by 
a mean-square displacement of the form ((Ax) 2 ') ~ t, and for the case of anomalous diffusion or 

subdiffusion, characterized by a mean-square displacement of the form ^(Aa;) 2 ^ ~ t" 1 with < 7 < f . 
In the context of classical diffusion, we obtain an expression for the mean first passage time and show 
that this quantity changes when the direction of the sawtooth is reversed or, equivalently, when the 
reflecting and absorbing boundaries are exchanged. We discuss at which numbers of 'teeth' TV (or 
number of DNA nucleotides) and at which heights of the sawtooth potential this difference becomes 
significant. For large N, it is well known that the mean first passage time scales as N 2 . In the 
context of subdiffusion, the mean first passage time does not exist. Therefore we obtain instead 
the distribution of first passage times in the limit of long times. We show that the prefactor in 
the power relation for this distribution is simply the expression for the mean first passage time in 
classical diffusion. We also describe a hypothetical experiment to calculate the average of the first 
passage times for a fraction of passage events that each end within some time t* . We show that this 
average first passage time scales as N 2 ^ 1 in subdiffusion. 



I. INTRODUCTION 

Recent studies on the passage of DNA through narrow 
channels, apart from being very interesting in their own 
right, have been motivated by the exciting possibility of 
developing a practical technique to characterize and se- 
quence DNA. These single molecule experiments were pi- 
oneered by Kasianowicz, Brandin, Branton and Deamer 
using a protein called a-hemolysin as a pore embed- 
ded in a lipid membrane. When voltage is applied across 
the membrane, ion current can be detected. When the 
DNA chain is inside and blocking the channel, the cur- 
rent is suppressed. By measuring the blockage time one 
can gather information about the voltage driven dynam- 
ics of DNA translocation. A large amount of exciting 
data have been accumulated through such experiments 
[H El- [J 01 ■ Recently, another group of experiments was 
performed using solid state nano-pores for voltage-driven 
DNA translocation 

On the theoretical side, much effort was placed in un- 
derstanding the entropic barrier associated with translo- 
cation of a long polymer 0, Q . The electrostatic barrier 
associated with charged DNA penetration through the 
low dielectric constant membrane has been recently dis- 
cussed theoretically Another interesting question 
is about friction: whether friction on the small piece of 
DNA passing through the narrow channel is larger or 
smaller than the friction experienced by the large DNA 
coils outside the membrane. To this end, Lubensky and 
Nelson assuming the channel friction dominance, 
were able to account for the observed bimodal distribu- 
tion of the passage times for single-stranded DNA. In 
their model, each base preferentially tilts towards one 
end of the single stranded DNA chain (5' or 3'. See fig- 
ure Each peak in the distribution of passage times 



therefore corresponds to a particular end entering the 
channel first, the channel being just wide enough for a 
single strand to pass through. More recently, Kantor 
and Kardar 

El 

and then Storm et al Ml argued that 
for sufficiently long DNA there must be a crossover to 
the regime of domination of the out-of-channel friction, 
in which case translocation should become subdiffusivc 
in character (with displacement growing slower than i 1 / 2 
with time t). 

A new spin is added to the story by the experiment 
by Meller and coworkers These authors devised 

an experiment involving DNA having a string of iden- 
tical bases (adenine) in the single-stranded portion and 
a hairpin loop at one end (figure^ held in place by bond- 
ing of complementary bases. Like double-stranded DNA, 
the hairpin cannot enter the channel (a transmembrane 
pore of a-Hemolysin) . The hairpin therefore constrains 
the DNA to enter the pore with its single-stranded end, 
as well as preventing the entire DNA from crossing the 
membrane. In their experiment, the DNA, driven by an 
applied voltage, enters the pore with its single-stranded 
end. Thereafter, once current is blocked by the DNA, 
the voltage is either switched off, in which case the DNA 
diffuses freely (non-driven), or the sign of the voltage is 
flipped, in which case the DNA is pushed back. More- 
over, by making two DNA samples, with the hairpin loop 
at opposite ends, it is possible to observe DNA sliding 
away from the pore in two opposite directions along the 
DNA contour, and the observation suggests that DNA 
escapes in one direction faster than in the other. 

In the experiment [Lll ]. by measuring the so-called 
"survival probability" S(t), which is the probability that 
a DNA molecule will stay in the pore as a function of the 
waiting time, it was determined that the voltage-free dy- 
namics of the 3' threaded molecules is about two times 
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slower than the corresponding diffusion of 5' threaded 
molecules having the same sequence. Importantly, in 
both cases the DNA was threaded from the same side of 
the pore (called the cis-side of a-HL). To delineate the 
underlying mechanism responsible for the observed dy- 
namics, the authors of the work 14] performed all-atom 
molecular dynamics simulations, which independently 
confirmed the experimental results for driven DNA. The 
simulations also showed that the confinement of the DNA 
bases in the a-HL pore results in an even stronger (com- 
pared to a free DNA) tilt of the bases with respect to the 
DNA backbone towards the 5' end. 

Authors of the work 0] phenomenologically interpret 
their data by assigning two different diffusion constants 
for the two separate experiments in which the same DNA 
is placed in the channel in two possible orientations. This 
interpretation is justified by the fact that the interac- 
tions between the DNA bases and the pore are different 
in these two cases (perhaps via different barrier heights 
within the framework of a sawtooth potential landscape 
discussed below). 

There is a temptation to summarize the experimental 
findings of the work 0] in one sentence (although no one 
made this mistake, including ^ij): DNA diffuses in one 
direction faster than in the other. Indeed, the observed 
asymmetry of dynamics is consistent with the tilt of the 
nucleotides with respect to the main DNA chain. This 
asymmetry then seems easy to understand if the anal- 
ogy is made with petting a cat along or against the grain 
of its fur; the cat responds very differently in the two 
cases (presumably because it experiences very different 
friction). Another, possibly even more obvious, analogy 
would be carrying a Christmas tree top first or base first 
through a narrow door; one again encounters very differ- 
ent resistance in the two cases. The point is that such 
analogies and interpretations are only possible for the 
driven DNA motion, particularly for the system far into 
the non- linear regime (in terms of force- velocity relation) , 
whereas for the portion of experiments in the work j 1 31 ] 
involving freely diffusing DNA such analogies and inter- 
pretations would be wrong; it is not surprising then that 
the authors of the work 14] did not use such analogies 
and interpretation for the freely diffusing DNA. Indeed, 
for free diffusion, the friction coefficient (averaged over 
the scale well exceeding a single base) moving in one di- 
rection and in the opposite direction must be the same, 
as follows from the Onsager symmetry relation, and the 
assumption of asymmetric friction would be a grave mis- 
take. Although no one actually made this mistake, in- 
cluding |14| . it is worth emphasizing why an assumption 
of asymmetric friction would be a mistake. Indeed, if we 
only imagine that DNA (not driven by any applied volt- 
age!) diffuses in one direction faster than in the other, 
then we can easily build a perpetuum mobile (see figure 
[2J moving indefinitely long through time at the expense 
of thermal energy from the thermal bath, which is, of 
course, impossible. In other words, freely diffusing DNA, 
when it is already in the pore, in contrast to a (heavily 



driven!) Christmas tree through a door, must have the 
same friction coefficients when the DNA moves in either 
direction. 

What is nice is that the experimental findings and their 
interpretation in the work [l4| are in fact in perfect agree- 
ment with this thermodynamic analysis. In order to 
make this reconciliation very clear, we immediately re- 
fer to the symmetry analysis in figure ^ Notice that the 
pore itself is asymmetric (its crystallographic structure is 
known ^H), the DNA backbone is also asymmetric (from 
3' end to 5' end) , and the loopy end creates further asym- 
metry. This gives four possible orientations of the pore 
and the DNA with the loop: two possibilities arise from 
two different mutual orientations of the DNA backbone 
with respect to the pore (indicated by the numbers 1 and 
2 in figure^]), and for each of these two orientations there 
are two possibilities to place the blocking loop (indicated 
by the letters A and B in figure QJ. This symmetry anal- 
ysis, as shown in figured is reminiscent of the symmetry 
analysis in the paper [ll| , except we have no electric field, 
but instead have loops at the DNA ends. 

We can now say that in any one of the arrangements, 
from 1A, IB, 2A or 2B, the DNA must experience the 
same friction moving up or down the pore; friction going 
up equals friction going down. At the same time, the 
friction in configurations 1A or IB can be different from 
friction in configurations 2 A or 2B, and they are likely 
to be different. That is why the work 0] assigns two 
different diffusion constants to the two DNA-pore mutual 
configurations (1A and 2A). By contrast, the loop itself 
likely has no effect on the friction or diffusion coefficient, 
so we expect that the diffusion coefficient should be the 
same for configurations 1A and IB (same goes for 2A and 
2B). In other words, there should be two distinct diffusion 
coefficients, not four. We shall argue in this work that, 
nevertheless, there will be four different diffusion times 
corresponding to the four configurations in figure ^ 

To explain our approach, it is convenient to adopt a 
terminology in which, instead of considering diffusion of 
the DNA chain, we consider diffusion of the passage point 
along the DNA contour. Following Lubensky and Nel- 
son 111, we consider a simple model in which asymme- 
try is presented in the underlying potential landscape. 
For simplicity, we model it with a sawtooth profile. The 
two orientations of the asymmetric potential (relative to 
the boundary conditions) correspond to the two possible 
placements of the blocking loop for a given orientation 
between DNA and pore (e.g. 1A and IB in figure QJ. 

Like Lubensky and Nelson ^J, we focus on the first 
passage time, which is the time it takes for the ini- 
tially fully- 'plugged' DNA to completely 'unplug' from 
the pore. In other words, it is the time needed for the 
diffusing particle (or a random walker) to arrive for the 
first time on the open end of the DNA, or to one end of 
the (0, L) interval, provided that a reflecting boundary 
condition is imposed at the opposite end. 

We would like to emphasize the fundamental differ- 
ence between asymmetric diffusion, which is prohibited 
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by thermodynamics, and symmetric diffusion over the 
asymmetric potential landscape. It is well known, and we 
show it explicitly in appendix[E| that stationary diffusion 
remains symmetric despite the asymmetry of the under- 
lying potential landscape, thus making nonfunctional the 
perpetuum mobile design of figure 

In this paper, we compute the mean first passage times 
(MFPT) corresponding to cases 1A and IB in figure^^r 
to cases 2A and 2B). We consider the brownian motion of 
a particle diffusing classically in an asymmetric sawtooth 
potential U(x) and in the inverted or reversed version of 
the potential (figure^. This model neglects the entropic 
barrier (of order UbT In N) presented by the DNA coils 
on both sides of the pore but through the consider- 

ation of subdiffusion it does take into account the extra 
friction created by those coils 01 ■ From the results, we 
discuss when the difference between the two times is sig- 
nificant. (Note that since we know very little about the 
details of the interactions between the DNA bases and 
the pore, we cannot determine if case A in figure [3] cor- 
responds to case 1A in figure ^ an( i case B in figure |3| 
corresponds to case IB in figure ^ or if it is the other 
way around.) 

Since DNA translocation is ultimately not classical dif- 
fusion, but rather subdiffusion 0, 0] , we consider also 
the first passage times for the subdiffusion in the presence 
of an asymmetric potential. In general, the first passage 
time for subdiffusion was recently a matter of consider- 
able interest and dispute in the literature 

ehgi it 

is now understood (Hi, I20L l2l| that the mean first pas- 
sage time diverges for subdiffusion, because a subdiffus- 
ing walker tends to remain too long on the place that 
it once reached. Accordingly, we look at the probabil- 
ity distribution for the first passage times (DFPT), and 
concentrate on its tail at long times. We found that this 
tail is very different for the two potentials, and the dif- 
ference turns out to be expressed through corresponding 
mean first passage times for classical diffusion. With this 
knowledge, we construct an average first passage time 
from a subset of passage events and show that this aver- 
age scales as N 2 / 1 . The result also exhibits the asymme- 
try between cases 1A and IB (or 2A and 2B) just as in 
the case of classical diffusion. 

From our discussion we make the prediction that the 
observed passage times for the four possible mutual ori- 
entations of the pore and the DNA will all be different. 



II. RESULTS 

A. Classical diffusion 

For treating classical or normal diffusion, one often 
starts with the Fokker-Planck (FPE) equation (in this 
context also frequently called Smoluchowsky equation) 



dP( x ,t) _ D d c _ 

dt dx 



^^■e^P^t) (1) 



giving the time-evolution of the probability density 
P(x,t). Here D is the usual diffusion constant and we 
have set ksT = 1. The FPE yields the Boltzmann dis- 
tribution for P in the steady-state, as well as giving the 
linear relation between the mean-squared displacement 
and time in the absence of external forces. 

In calculations involving the first passage time, it 
would be convenient to consider the equivalent problem 
of first passage to either x = L = Na or x = —L = —Na, 
where the potential U(x) for x > is as illustrated in fig- 
ure|31 while the potential for x < is U(x) for positive x 
reflected about the vertical axis. With this picture, the 
probability for the particle to still be 'alive' at time t, 
also called the survival probability S(t) (and measured 
in experiment 14]), is given by S(t) — j\P{x,t)dx. 
The distribution of first passage times F(t) is calculated 
from S(t) via F(t) = -^jj^. This gives the following 
expression for the mean first passage time t(xq) |25j : 



r(x ) 



tF(t)dt 



S{t)dt 



OO pL 



P(x, t)dxdt 



(2) 



where xq is the initial position of the particle, P(x, 0) = 

S(x - X ). 

It can be shown that r(xn) satisfies an ordinary dif- 
ferential equation [2^, 0, El (derived in appendix [5J| • 
The solution of this differential equation for a sawtooth 
potential U(x) is outlined in appendix IH1 For the parti- 
cle initially located at the origin (x$ — 0), the mean first 
passage time to reach x = L = Na is given by, 



ta 



L 2 ah 
a— — p- 



2D ' D 
for the potential in figure [3J\, and 

L 2 aL 

tr = a 1- p — 

2D D 



(3) 



(4) 



for the potential in figure [3)3. Here, we have defined the 
coefficients 



a 







^mh{U /2) 



U /2 
sinh (Uq) - 



U 



JT 2 



(5) 
(6) 



Expression (@J can be obtained from © by flipping the 
sign of Uq. 

From the results © and Q, it is clear that ta < tb- 
Physically, the inequality ta < tb may be obvious for 
the case N = 1 in figure |3 in which a particle has to 
surmount a single barrier in order to get to x = L in 
case B, while there is no barrier in case A. In general, 
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for a given N, the potential in figure involves N — 
1 barriers, while the potential in figure |3j3 involves N 
barriers. In fact, it is easy to show that in the limit 
Uq 3> 1 (Ua/ksT ^> 1 in more conventional units) we 
have ta(N + 1) — tb(N), where the arguments indicate 
the number of teeth in the sawtooth potentials. 

For the long DNA, when L > a or iV > 1, the leading 
terms in both ta © and tb © are proportional to A 2 , 
as one would expect for diffusion times. To this leading 
order, first passage times ta and tb obey the symmetry 
in diffusion and are the same. It is in the subleading 
terms (proportional to L) that the two times differ. Let 
us stress that the difference between ta and tb , which is 
of order of \/N in a relative sense, is entirely due to the 
boundary conditions and the situation at the ends of the 
diffusion region. 



B. Anomalous diffusion 



Anomalous diffusion is characterized by the occurrence 
of a mean square displacement of the form ((AxY) ~ V , 
where < 7 < 1 in subdiffusion; traditionally [22j. this 
is written in the form 



<(A,) 2 ) 



2D-, 



r(i + 7) 



(7) 



where _D 7 is a generalized diffusion constant and T(x) is 
the gamma function. For 7 = 1 one recovers the usual re- 
sult for classical diffusion. It can be shown that this 
form for the mean square displacement can be obtained 
from a generalized version of equation © called the frac- 
tional Fokker-Planck equation (FFPE). This equation is 
described in appendix 1X1 

Although up to this point we have ignored the inter- 
actions of the DNA bases outside the pore, it seems rea- 
sonable to speculate that their effect is to slow down the 
translocation. Thus one might be able take these interac- 
tions into account phenomenologically by positing a value 
of 7 corresponding to the subdiffusive domain < 7 < 1. 
(Reference 0] lists possible sources of waiting time dis- 
tributions leading to anomalous diffusion ). 

It is shown below and in references 

0G1HII1 that 
the MFPT does not exist for subdiffusion. This leads us 
to consider the probability distributions themselves. The 
method of Laplace transforms can be used to solve for the 
transform of the survival probability [l6l Hsj but one is 
left with the very difficult task of obtaining the inverse 
transform, even for the case of a sawtooth with N = 1. 
However, it is shown in appendix \B\ that the long-time 
limit of the survival probability and the first passage time 
distribution scales as some power of t and that they are 
simply related to the expression for the MFPT in the 
context of classical diffusion as follows 



S{t) 



T(l- 7 ) 



(8) 



F(t) = 



dS(t) 
dt 



(9) 



Here r 7 is the same expression as the MFPT in classi- 
cal diffusion © (or Q), but containing a generalized 
diffusion constant. The long-time limit is reached when 
£T > t 7 ~ L 2 /Dj ~ A 2 . The relationships © and © 
ultimately arise from the almost identical expressions for 
the solution P(x,t) in classical diffusion and in subdif- 
fusion. The two solutions differ only in the time depen- 
dence, which is an exponential for classical diffusion. 

Because the expectation time of a first passage is infi- 
nite, any meaningful experiment, real or computational 
alike, must be based on some protocol rendering the ob- 
servation time finite. We argue that in essence such a 
protocol is always reduced to discarding the events which 
fail to come to completion within some specified time t* ; 
in other words, only those passage events that each com- 
plete within some time t* are counted. The rest of the 
events that do not end by time t* are terminated and dis- 
carded. The conditional probability distribution of the 
first passage events that get counted under such a proto- 
col is then given by 



F(t) 



l-f t 7F(t)dt 



(10) 



For such an experiment, there exists a perfectly defined 
and finite average first passage time. This conditional 
average, for large t*, is 



Jo tF(t)dt 
l-f t 7F(t)dt 



it" 



'7(l- 7 )r(l- 7 ) 



1 



t*TT(l- 7 ) 
7 t* 1 ~ 7 

(i- 7 )r(i- 7 ) 



rT(i- 7 ) 



(ii) 



(12) 



So far, the time t* should be long enough, but otherwise 
arbitrary. Now we argue that the time t* must be chosen 
such that roughly about half of the passage events at a 
given A get discarded. This requirement seems reason- 
able, for if one discards a much smaller fraction t* be- 
comes too large and the measurements get inefficiently 
slow; if one discards a much larger fraction t* becomes 
too small and the tail of the distribution does not get 
sampled properly. Thus, assuming half of the events dis- 
carded, t* becomes of order (r 7 ) 1 / 7 , just at the boundary 
of the validity of the asymptotics. Substituting this into 
(|llf> . one obtains a scaling of (tj) 1 ^ 1 ~ TV 2 / 7 for the av- 
erage first passage time. Of course, this scaling is not un- 
expected for subdiffusion with an average displacement 
going like t 1 ! 2 . Furthermore, due to the appearance of 
the classical diffusion times ta and tb (which take the 
place of ^ depending on the potential) in the average 
first passage time we just defined, the asymmetry of the 
first passage time is once again present in this case. 
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III. DISCUSSION 



we find an average that scales as (r 7 ) 1/7 ~ N 2 /~t '. 



The ratio of the MFPTs in classical diffusion, expres- 
sions (01 and is plotted in figure 01 for a few realistic 
values of N and Uq. We see that for Uq equal to a few 
ksT, the difference becomes small (~ 10%) for N > 10. 
For N = 50, corresponding to the length of ssDNA used 
in the experiments by Meller and coworkers and 
for U a /k B T ~ 10, the fractional difference in MFPTs is 
about 4% . 

We emphasize again that one cannot use the results of 
the comparison between these two times (cases 1A and 
IB in figure^) and apply it to the experimental results in 
[l4| (cases 1A and 2A in figure^). Due to the asymmetry 
of the pore, the 3' and the 5' threading of DNA through 
one end of the pore (the so-called cis side) cannot be 
readily reduced to cases 1A and IB in figure 

Having established the difference in average first pas- 
sage times for the two asymmetric potentials, let us now 
turn to the scaling of the first passage times with TV. For 
large N, the scaling result N 2 found earlier is well known 
for classical diffusion or Brownian dynamics. However, 
this is in conflict with the equilibration time of a poly- 
mer with N monomers in the absence of a pore and 
membrane, which already scales with N to some power 
larg er than 2 for Rouse dynamics of self-avoiding chains 
|13| . This suggests that a correct description of polymer 
translocation should be made in the context of subdiffu- 
sion, where the scaling of the average first passage time is 
to power I/7 > 1 of the classical result, although we do 
not give a prediction for the value of 7 itself because the 
interactions involving the DNA/polymer located outside 
the pore were not treated explicitly. The scaling TV 2 / 7 is 
also not surprising if one takes the relation ((Ax) 2 ) ~ t 1 
and puts Ax ~ N, but it does not rule out the argument 
made above regarding t* , only that it gives a reasonable 
and somewhat expected answer. Moreover, using similar 
arguments, the result N 2 ^ 1 is consistent with numerical 
simulations made by Chuang, Kantor and Kardar |13| for 
diffusive dynamics of self- avoiding chains in two dimen- 
sions. They found that the average of the first passage 
time scales as r ~ N 2 - 5 = iV 1+2ly , where v = 3/4 in 
two dimensions. They also argued, assuming that the 
translocation coordinate goes like (Aa: 2 (i)) ~ i 7 at short 
times, that 7 = 2/(1 + 2v). Eliminating 1/, their formulas 
imply that r ~ N 2 ^ . 

To summarize, based on our calculations for the mean 
first passage times in asymmetric sawtooth potentials 
and experiments by Meller and coworkers [T^j , we expect 
that the average first passage times for the four cases in- 
dicated in figure^are all different. The expression for the 
tail of the first passage time distribution in subdiffusion 
is of the form jTy/T(l — 7)i 1+7 , where r 7 is the formula 
for the mean first passage time in classical diffusion. Be- 
cause the power of t in the distribution is less than two, 
the mean first passage time diverges. By constructing 
an average from the first passage times less than time t* 
such that approximately half of the passages get rejected, 
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APPENDIX A: FRACTIONAL 
FOKKER-PLANCK EQUATION 

A generalization of the FPE describing anomalous dif- 
fusion is given by the fractional FPE |22j 



dP(x,t) 
dt 



oD 



1-7 i 



TT 



P 



Equivalently, 



t-'P(x,0) 
oD t 'P{x,t) — — = L 



FP 



P 



r(i- 7 ) 

where the Fokker-Planck operator is defined as 
FP 7 ftr dx 



(Al) 



(A2) 



(A3) 



Here Z? 7 is a generalized diffusion coefficient and U (x) 
is an external potential. We have also set fcsT = 1 and 
the Einstein relation is implicit. The Riemann-Liouville 
fractional operator is defined through 



1 d 



r( 7 ) dt J (t - 1') 1 - 7 



(A4) 



One can easily check that the FFPE reduces to the FPE 
or diffusion equation for 7 = 1. 

Given the initial distribution P(x,0) = 5(x — Xq), the 
solution to equation i|Alfl is given by the bilinear expan- 
sion El 



P(x,t;x ,0) 



e U(x„)/2-U(x)/2 x 



Y.M^nME^-K^) . (A5) 



The functions <f> n (x) = e^ u ^/ 2 ^ n {x) and T„(i) = 
E~ f {— A„i 7 ) appear in the separation of variables ansatz 
W n (x,t) = 4> n (x)T n (t). The product function W n (x.t) 
satisfies the FFPE. Note that the coordinate dependence 
comes through the eigenfunctions ip n (x) or <fi n (x), which 
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are the same as for regular diffusion, satisfying the (eigen- 
value) equations 



L-pp^n(x) = -X n 4>„(x) , 



(A6) 



Taking the inverse transform, one obtains the long-time 
behaviour of the Mittag-Leffler function 



4p m,t "i(i) 



-X n ip n (x) , 



(A7) 



£ 7 (-At 7 ) 



(-1) 



m+l 



' r(i — 7m) 

m=l y ' ' 



(A13) 



L ^er,nitian _ & V '(x) /2 L^-V '(*) /2 _ ( Ag ) 

However, as to the time dependence, which for classical 
diffusion is described by exponentials (e _A ' lt ), for subd- 
iffusion it must satisfy the equation 



dT n {t) 
dt 



-A„ oD t T n (£) 



(A9) 



One can check that the following series definition of the 
Mittag-Leffler function E 1 (z) satisfies equation (|A9|) 



r(i + 7 m) 



(A10) 



This function is a natural extension of the exponential 
function, to which it degenerates for 7 = 1. 

By taking the Laplace transform of both sides of equa- 
tion 1|A9J| . one obtains an alternative definition of the 
Mittag-Leffler function 



£{£ 7 (-At 7 )} = (s + Xs 1 - 1 ) 1 



(All) 



(The subscript in the constant A has been dropped.) The 
long-time limit of the Mittag-Leffler function corresponds 
to the small s limit of the Laplace transform. Expanding 
(|A11|) in a series for small s, 



For e=l—7 = 0, the Laplace transform IjAllfl be- 
comes (s + A) -1 , the inverse transform of which is an 
exponential. For e close to 0, we expect a long time 
interval in which the behaviour of the Mittag-Leffler 
function E 7 {— Xt 1 ) behaves like an exponential; at much 
longer times the behaviour changes to a power law. The 
crossover is expected to happen when e 



-At c 



t, or 



r(e)At c ' 
at about t c ~ j ln(-)- 

APPENDIX H: DIFFERENTIAL EQUATION 
SATISFIED BY THE MEAN FIRST PASSAGE 
TIME 



Recall that the MFPT can be calculated from (equa- 
tion ©) 



r(x ) 



00 pL 



P{x, t)dxdt 



(Bl) 



C{E 7 (-X^)} 



1 



As 1 " 7 



1 "T + T 1 - 



rn — 1 



,7771-1 

A m 



(A12) 



To derive an ordinary differential equation satisfied by 
t(xo), apply the operator e u ' x °^Lpp ^ e - u (, x o) to equa- 
tion (|B1(I and use the eigenfunction expansion solution 



e U{xo)L FF, Xo e ~ U{x °Mxo) = 



J-L 
L 



^/^^/^(-^M^M^E^-x^dxdt 

t-tP(x,Q) 



71=0 



J-L 



Lyy>P( x i t)dxdt 



J-L 



D?P(x,t) 



r(i- 7 ) 



dxdt 



In the last two steps, the eigenvalue equations and the 
second version of the FFPE (equation (|A2|i 1 was used. 

Using the initial condition P(x, 0) = S(x — xq) and the 
definition of the fractional operator, after some algebra 



one obtains 

e^Lpp^e-^VOco) = 

t-00 [r(2- 7 ) r(i- 7 )7 (t-t'y \ y ' 



7 



D e u iX0 ) A_ e - 

t l--i 



ox 



— lim 

t — >OC 



r(2- 7 ) r(i- 7 )y 



-eft' 



(B3) 



For 7 = 1, corresponding to classical diffusion, the sur- 
vival probability S(t') decays exponentially and the term 
with the integral goes to zero, yielding the familiar result 
of -1 for the right-hand-side [H HI H^. For 7 < 1, 
S(t') goes like (£')~ 7 (see (|D2j|) and the term with the 
integral goes like i 1-27 . The right-hand-side diverges, 
which hints at the non-existence of the MFPT for subd- 
iffusion [H El El HI EH- 



APPENDIX C: SOLUTION FOR THE MEAN 
FIRST PASSAGE TIME IN A SAWTOOTH 
POTENTIAL 

From the previous section, the differential equation 
satisfied by the MFPT in the context of classical diffusion 
is (temporarily putting back fcsT) 



Dc U(x)/k B T d c -U( x )/k B T dT ( x ) = _j 



dx 



dx 



(CI) 



We solve for t[x) in this equation for a sawtooth po- 
tential (case A, figure [2J subject to the boundary con- 
ditions dr/dx(0) = and t(L) = 0, and the continuity 



of r and e v l kBT dr j dx in (0,L). In what follows we let 



a va/D 



oU /k B T 



The solution, for x between (m — l)a and ma where m 
is an integer between 1 and N (inclusive), is given by 



r(x) = A m e- VX ' D — — + B. 



(C2) 



The coefficients A and B are given by 



A m = r- 1 ^ [(£-l)(m-l)-l] 

D (a , D 
B m = — + NI-- 1-1/0- 

D{i - l) 2 (N - m)(N + m - 1) 



(C3) 



The MFPT for a particle initially located at x = is 
given by r(0) = Ai + B±. To obtain the solution for case 
B in figure^ we may flip £ (— > l/£) and the sign of i> in 
the expressions above. 



APPENDIX D: RELATIONSHIP BETWEEN THE 
DFPT IN ANOMALOUS DIFFUSION AND THE 
MFPT IN CLASSICAL DIFFUSION 

Since the MFPT does not exist for subdiffusion, one 
would want to calculate the distributions instead. In the 
long time limit, using JXHJ and JUI!, 



lim P 7 (x,t ;a:o ,0)^^e^/ 2 -^/V»(^M rn \. - 
t-00 l(l-7)A n n 



(Dl) 



lim 5 7 (t) 



E 



L / v uy r(i- 7 )A„n 



(D2) 



r 



(D3) 



Again, these results indicate that the MFPT diverges 
for 7 < 1. It is also interesting to note that all the 
eigenfunctions ip n (x), not just the ground state, enter in 
the expressions. 

To make sense of the expression multiplying r ^_^ n+1 
in (|U3|) write down the corresponding solution for classi- 
cal diffusion under the same potential and the same value 
for the diffusion coefficient 

P(x,t;X ,0) = e U(*a)/2-U( X )/2 x 



x E tpn{x)ip n (xo)exp(-\ n t) (D4) 

(Note exponential instead of Mittag-Leffler function). 
The survival probability is given by 

S(t) = ^(j^'^^M*)**) X 

x e^/ 2 V>„(xo)exp(-A^) (D5) 
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While the MFPT is given by 

t(x ) = / S(t)dt 
Jo 



X e ^o)/^ n(a;o) J_ 



(D6) 



which is identical to the coefficient of 7 /r(l — "f)t 1+1 in 



aries. Let the potential f7(x) satisfy ?7(0) = f/(£) = 0, 
but is otherwise arbitrary. The classical diffusion equa- 
tion is given by 



dc 

at 



dJ_ 

dx 



(El) 



where J = — De~ u ^ -§^& u ^ x 'c (see equation (Q). In the 
steady state, §§ = 0, which implies that J is spatially 
uniform. Integrating Je u ^ = —D-S^e u ^ x 'c{x) and uti- 
lizing the boundary conditions, one obtains 



APPENDIX E: EFFECTIVE DIFFUSION 
CONSTANT IN THE STEADY STATE 



J = 



D 



ft e u Mdx 



(c(O)-c(L)) 



(E2) 



In this section we determine the steady state current 
J given fixed concentrations c(0) and c{L) at the bound- 



This expression is identical with Fick's law with an effec- 
tive diffusion constant of nT — . 

Jo 
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FIG. 1: (Color online) The four possible relative orientations 
of DNA ( "key" ) and pore ( "keyhole"). In a similar figure from 
Lubensky and Nelson (reference figure 7), the single- 

stranded DNA has no loop. Instead, the four cases were due 
to the various relative orientations of DNA, pore and an ap- 
plied electric field. Arrows in our figure show the direction 
from the 5' to the 3' end in the DNA. Inset in the middle 
shows schematically the tilted bases. The analysis in our work 
compares the passage times for case 1A with case IB (or 2A 
with 2B) in which the relative orientation between the DNA 
bases and pore is identical. In contrast, the experiments in 
study cases 1A and 2A, where the DNA enters the pore 
from the same (cis) side. 



FIG. 2: (Color online) This arrangement of DNA and pores 
would have acted as a perpetuum mobile if the stationary dif- 
fusion coefficient was asymmetric. This shows that it cannot 
be asymmetric, the symmetry being a requirement of thermo- 
dynamics. 
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FIG. 3: (Color online) The sawtooth potential of N teeth with 
a reflecting boundary at x = and an absorbing boundary at 
x = L = Na, illustrated for two different directions of asym- 
metry in the sawtooth. The reflecting boundary corresponds 
to the inability of the DNA hairpin to pass through the pore. 
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FIG. 4: (Color online) The ratio ta/tb plotted against the 
dimensionless drift or barrier height va/D = Uo/kBT for N = 
1 (bottommost curve), N = 2, N = 5 and N = 10. 



